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The ion-sputtering induced intermixing is studied by Monte-Carlo TRIM, molecular dynamics 
(MD) simulations, and Auger electron spectroscopy depth profiling (AES-DP) analysis in Pt/Ti/Si 
substrate (Pt/Ti) and Ta/Ti/Pt/Si substrate (Ti/Pt) multilayers. Experimental evidence is found 
for the asymmetry of intermixing in Pt/Ti, and in Ti/Pt. In Ti/Pt we get a much weaker interdif- 
fusion than in Pt/Ti. The unexpected enhancement of the interdiffusion of the Pt atoms into the Ti 
substrate has also been demonstrated by simulations. We are able to capture the essential features 
of intermixing using TRIM and MD simulations for ion-beam sputtering and get reasonable values 
for interface broadening which can be compared with the experimental measurements. However, 
the origin of the asymmetry remains poorly understood yet. 
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I. INTRODUCTION 

Bombardment of surfaces by energetic particles of- 
ten leads to atomic inter-layer mixing [1-6]. The ion- 
induced inter-layer atomic transport through interfaces 
(ion-intermixing) has been the subject of numerous stud- 
ies [2,7,8]. However, recent results indicate that our 
knowledge on the mechanism of ion-mixing might be in- 
complete [9]. 

Ion-sputtering of surfaces by energetic particles has 
widely been used for depth profiling of interfacial struc- 
tures [10]. Recent advances in the use of ion-sputtering 
has been attracted considerable theoretical and experi- 
mental attention [11]. In particular, ion-sputtering has 
been simulated by molecular dynamics using sequential 
bombardment of metallic surfaces [12-15]. In order to 
get more insight into processes which govern interdiffu- 
sion during repeated ion-bombardments Auger electron 
spectroscopy (AES) depth profiling analysis (see e.g. the 
review refs. of [16]) has been used for the study of broad- 
ening of interfaces during ion-sputtering in combination 
with transmission electron microscopy (TEM) [1,17]. 

There is now also a widespread interest in understand- 
ing the asymmetry and the anomalously wide range of in- 
termixing in various diffusion couples [18,19]. It has also 
been reported that interdiffusion is not driven by bulk 
diffusion parameters nor by thermodynamic forces (such 
as heats of alloying) [21,18]. Computer simulations have 
also revealed that mass-anisotropy of the bilayer governs 
ion-bombardment induced interdiffusion at the interface 
[9] and greatly influences surface morphology develop- 
ment during ion-sputtering [24,25]. We also pointed out 
that interdiffusion takes place via ballistic jumps (ballis- 
tic mixing) in various diffusion couples under the effect 
of ion-bombardment [21]. 

We demonstrate in the present communication that 
using atomistic simulations we are able to reproduce the 
experimental interfacial broadening and the asymmetry 



of intermixing obtained by ion-sputtering. Furthermore, 
we would like to model and explain the enhancement of 
intermixing as a function of the anisotropy of the inter- 
faces. We reproduce and partly explain the experimen- 
tally found asymmetry of intermixing by TRIM and MD 
simulations. 



II. THE SETUP OF THE MEASUREMENTS AND 
THE SIMULATIONS 

A. The experimental setup 

The experimental setup is the following: According to 
the crossectional TEM (XTEM) results the thickness of 
the layers in samples are: Pt 13 nm/Ti 11 nm/Si sub- 
strate (denoted throughout the paper as Pt/Ti), and Ta 
21 nm (cap layer to prevent oxidation of Ti) /Ti 11 nm/Pt 
12 nm/Si substrate (denoted as Ti/Pt). The XTEM im- 
ages are shown in Fig 1. For the sake of simplicity we 
consider our multilayer samples as bilayers and we study 
the atomic transport processes at the Ti/Pt and Pt/Ti 
interfaces. Both samples have been AES depth profiled 
by applying various sputtering conditions. The sample 
has been rotated during sputtering. In the following we 
will outline results of 500 eV Ar + ion bombardment at 
an angle of incidence of 10° (with respect to the sur- 
face). The atomic concentrations of Pt, Ti, Ta and Si 
were calculated by the relative sensitivity method taking 
the pure material's values from the spectra. The oxygen 
atomic concentration has been calculated by normaliz- 
ing the measured oxygen Auger peak-to-peak amplitude 
to TiC>2 [36]. The depth scale was determined by as- 
suming that the sputtering yield (Yi) in the mixed layer 
is the weighted sum of the elemental sputtering yields 
*V .Y,V)). 

The broadening of the interface is frequently charac- 
terizied by the depth resolution. The depth resolution 
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is defined as the distance of points on the depth profile 
exhibiting 84 % and 16 % concentrations. This definiton 
has been introduced for the "normal" cases, when either 
the ion mixing or the roughening can be described by a 
Gaussian convolution resulting in an erf fuction transi- 
tion in the depth profile. The same definition used for 
other cases as well, however. 

If the transition does differ from the erf function (e.g. 
when the mobility of one of the components of a diffusion 
couple is much higher than that of the other's) one might 
give also the distances between points of 84 % and 50 %, 
and 50 % and 16 %. The ratio of these distances gives 
us the asymmerty of intermixing (shown in Table 1). 

B. The computational methods 

Dynamic TRIM (Transport of Ion in Matter) TRIDYN 
simulation has been applied to model the ion sputter- 
removal process (see e.g. ref. [26]). The input parameters 
of the code used (TRIDYN [27]), characterizing the pure 
material, are: atomic number and density, (pn = 56.8 
at/nm 3 , ppt = 66.2 at/nm 3 ), bulk binding energies (0 for 
both materials), and the surface binding energies (SBE) 
of Ti (4.88 eV) and Pt (5.86 eV). The simulation provides 
the atomic concentration along the depth after a given 
dose of bombarding ions. Having this distribution the 
Auger intensities of elements present can be calculated 
using standard equations and parameters like inelastic 
mean free path (IMFP), backscattering factors, primary 
current etc (the software provided by S.Tougaard is used 
to calculate the IMFP based on TPP-2M). Repeating this 
procedure (simulation of concentration distribution after 
additional ion bombardment and calculation of the cor- 
responding Auger intensities) until the layer sputtered 
away results in a simulated depth profile, which might 
be compared with the measured one. (The complete de- 
scription of this calculation is given in ref. [28]). As will 
be seen later the code describes the essence of the experi- 
ments, but fails to provide a quantitative agreement. It is 
not surprizing since we know that this code accounts nei- 
ther for roughening nor for the recovery processes taking 
place after the collisional cascade. 

In order to get more insight to the mechanism of inter- 
diffusion classical molecular dynamics simulations have 
also been used to simulate the ion-solid interaction (using 
the PARC AS code [29]). Here we only shortly summarize 
the most important aspects. A variable timestep and the 
Berendsen temperature control is used to maintain the 
thermal equilibrium of the entire system. [30]. The bot- 
tom layers are held fixed in order to avoid the rotation of 
the cell. Periodic boundary conditions are imposed lat- 
crarily and a free surface is left for the ion-impacts. The 
temperature of the atoms in the outermost layers was 
softly scaled towards the desired temperature to provide 
temperature control and ensure that the pressure waves 
emanating from cascades were damped at the borders. 



TABLE I. The experimental and TRIDYN results for 
depth resolution and for the asymmetry of intermixing for 
Ti/Pt and Pt/Ti samples. 





exp 


TRIDYN 


asym. (exp) 


asym. (TRYDIN) 


Ti/Pt 


2.0 


3.3 


0.9 


0.9 


Pt/Ti 


7.0 


4.3 


2.3 


2.0 



The exp and TRIDYN denote our measured and calculated 
depth resolutions (nm) . In the 3rd and 4th columns the exper- 
imental and the calculated asymmetry of mixing are given. 



The lateral sides of the cell are used as heat sink (heat 
bath) to maintain the thermal equilibrium of the entire 
system [30]. The detailed description of other technical 
aspects of the MD simulations are given in [29,30] and 
details specific to the current system in recent commu- 
nications [9,21,24,25]. Atomic collisions have been sim- 
ulated in a standard way given in refs. [29,33]. Recoils 
have been initialized by giving the kinetic energy of the 
incoming ion to a lattice atom which is nearby the impact 
position. 

FIG. 1. The crossectional TEM (XTEM) images of the as 
received samples. Upper panel: Ta/Ti/Pt/Si, the thickness of 
the Ta, Ti and Pt layers are 21, 11 and 12 nm, respectively. 
Lower panel: Pt/Ti/Si, the thickness of the Pt and Ti layers 
are 13 and 10 nm, respectively. 



We irradiate the bilayers Pt/Ti and Ti/Pt with 0.5 
keV Ar + ions repeatedly with a time interval of fO-20 
ps between each of the ion-impacts at 300 K which we 
find sufficiently long time for the termination of intcrdif- 
fusion, such as sputtering induced intermixing (ion-beam 
mixing) [21]. The initial velocity direction of the impact- 
ing ion was 10 degrees with respect to the surface of the 
crystal (grazing angle of incidence) to avoid channeling 
directions and to simulate the conditions applied during 
ion-sputtering. We randomly varied the impact position 
and the azimuth angle <f) (the direction of the ion-beam). 
In order to approach the real sputtering limit a large 
number of ion irradiation are employed using automa- 
tized simulations conducted subsequently together with 
analyzing the history files (movie files) in each irradiation 
steps. In this article we present results up to 200 ion ir- 
radiation which we find suitable for comparing with low 
to medium fluence experiments. 200 ions are randomly 
distributed over a 2.0 x 2.0 nm 2 area which corresponds 
to <~ 5 x 10 15 ion/cm 2 ion fluence and to the removal of 
few MLs. 

The size of the simulation cell is 11.0 x 11.0 x 9.0 
nm 3 including 57000 atoms (with 9 monolayers (ML) 
film/substrate). At the interface (111) of the fee crys- 
tal is parallel to (0001) of the hep crystal and the close 
packed directions are parallel. The intcrfacial system is 
a heterophase bicrystal and a composite object of two 



2 



different crystals with different symmetry is created as 
follows: the hep Ti is put by hand on the (111) Pt bulk 
(and vice versa) and various structures are probed and 
are put together randomly. Finally the one which has 
the smallest misfit strain prior to the relaxation run is se- 
lected. The difference between the width of the overlayer 
and the bulk does not exceed 0.2 — 0.3 nm. The remain- 
ing misfit is properly minimized below ~ 6% during the 
relaxation process so that the Ti and Pt layers keep their 
original crystal structure and we get an atomically sharp 
interface. During the relaxation (equilibration) process 
the temperature is softly scaled down to zero. Accord- 
ing to our practice we find that during the temperature 
scaling down the structure becomes sufficiently relaxed 
therefore no further check of the structure has been done. 
Then the careful heating up of the system to 300 K has 
been carried out. The systems were free from any serious 
built-in strain and the lattice mismatch is minimized to 
the lowest possible level. The film and the substrate are 
~ 2.0 and ~ 6.8 nm thick, respectively. 

In order to reach the most efficient ion energy depo- 
sition at the interface, we initialize recoils placing the 
ion above the interface by 1 nm (and below the free sur- 
face in the 9 ML thick film) at grazing angle of inci- 
dence (10° to the surface) with 500 eV ion energy. In 
this way we can concentrate directly on the intermixing 
phenomenon avoiding many other processes occur at the 
surface (surface roughening, sputter erosion, ion-induced 
surface diffusion, cluster ejection, etc.) which weaken en- 
ergy deposition at the interface. Further simplification is 
that channeling recoils are left to leave the cell and in 
the next step these energetic and sputtered particles are 
deleted. 

We used a tight-binding many body potential, devel- 
oped by Cleri and Rosato (CR) on the basis of the second 
moment approximation to the density of states [34], to 
describe interatomic interactions. This type of a poten- 
tial gives a good description of lattice vacancies, including 
atomic migration properties and a reasonable description 
of solid surfaces and melting [34] . Since the present work 
is mostly associated with the elastic properties, melting 
behaviors, interface and migration energies, we believe 
the model used should be suitable for this study. The in- 
teratomic interactions are calculated up to the 2nd near- 
est neighbors and a cutoff is imposed out of this region. 
This amounts to the maximum interatomic distance of 
~ 0.6 nm. For the crosspotential of Ti and Pt we employ 
an interpolation scheme [9,24,33] between the respective 
elements. The CR elemental potentials and the interpo- 
lation scheme for heteronuclear interactions have widely 
been used for MD simulations [9,20,31,32]. The Ti-Pt 
interatomic crosspotential of the Cleri-Rosato [34] type 
is fitted to the experimental heat of mixing of the corre- 
sponding alloy system [21,25]. The scaling factor ro (the 
heteronuclear first neighbor distance) is calculated as the 
average of the elemental first neighbor distances. 

The computer animations can be seen in our web page 
[35]. Further details are given in [29] and details specific 



to the current system in recent communications [9,21]. 

III. RESULTS 
A. Experimental results 

Two typical AES depth profiles recorded on samples 
of Pt/Ti and Ti/Pt at eV ion energy, are shown, re- 
spectively in Fig. 2. It is clear without any detailed 
evaluation that the Pt/Ti and Ti/Pt transitions are very 
different. While the Ti/Pt transition is a "normal" one, 
in the case of Pt/Ti an unusual deep penetration of Pt 
to the Ti phase is observed. 

According to the experiment a relatively weak inter- 
mixing is found in Ti/Pt (<r « 2.0 nm) while an unusually 
high interdiffusion occurs in the Pt/Ti bilayer (<r « 7.0 
nm) . Moreover we observe a long-range tail in upper Fig 
2 for Pt in Pt/Ti while no such tail is found for Ti/Pt for 
Pt (neither for Ti, see lower Fig 2). The depth profiles 
also show oxygen presence as well. Does the presence of 
O influence the result? The answer is no which is detailed 
bellow. 

We have carried out several experiments, using differ- 
ent parts of the sample and applying various sputtering 
conditions. In all cases we recognized oxygen. The oxy- 
gen concentration varied considerably. The oxygen at the 
metal/substrate interface is due to the native oxide on the 
Si substrate. Since the bulk level of oxygen slightly cor- 
related with the ion current intensity (the larger the ion 
current intensity the lower the oxygen AES signal) part 
of the oxygen is contamination occurring during the AES 
depth profiling process. On the other hand the interface 
broadening does not correlate with the concentration of 
the oxygen in the Ti/Pt interface. Thus we conclude 
that the atomic transport is not affected by the presence 
of the slight oxygen contaminant. 

XTEM gives reasonable good estimate of the broad- 
ening of the as received sample, which was found to be 
about ~ 1,5 nm (intermixing would result in a gradual 
change of the contrast; since it was missing we estimate 
that the initial intermixing is less than ~ 1, 5 nm). Ion- 
sputtering might cause surface roughening as well as ion 
mixing. In case of our experimental conditions (rotated 
sample, grazing angle of incidence, low relative sputtering 
yield) the ion bombardment induced roughening is ex- 
pected to be weak, but cannot be ruled out; thus part of 
the measured total broadening might be due to roughen- 
ing. 

Unfortunately, using MD simulations it also hard to 
account for interface roughening, since the wavelength of 
interface roughening often exceeds the lateral size of the 
simulation cell. Therefore, we account only for intermix- 
ing, and then we estimate the magnitude of roughening 
as a difference of the measured broadening a and the 
simulated one. 
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FIG. 2. The concentration depth profile as a function of the 
removed layer thickness (nm) obtained by AES depth profiling 
analysis using ion-sputtering at 500 eV ion energy (10° with 
respect to the surface) in P t/Ti (upper Fig: 2a) and Ti/Pt 
(lower Fig: 2b). 



The measured depth resolution can be considered as 
the broadening of the interface a in those cases when the 
inelastic mean free path of the Auger signal electrons is 
much smaller than the depth resolution. This condition 
holds f or the case of Ti and Pt. Moreover, a is defined as 
a = V a ro + °mix ' where a ro and a mix are the roughen- 
ing and intermixing components [16]. Unfortunately we 
can not extract the components even in the as received 
samples, hence the cr w 1,5 nm measured by XTEM in 
the parent samples might include components both from 
interface roughening and intermixing. Also, we can say 
not too much about the ratio of oy to a m ix after ion- 
sputtering. We measure an overall value for broadening. 
We employ TRIM and molecular dynamics simulations 
to account for the contribution of intermixing a m i X to 
broadening a. 



B. Results obtained by MD and TRIDYN 

In Table 1 the experimental and TRIDYN results are 
summarized for the depth resolution and for the asym- 
metry of mixing. It is clear that the TRIDYN describes 
the essential features of the experiments. It predicts that 
broadening is wider when the Pt layer is on the top of 
the Ti, and that the penetration of the Pt into the Ti 
is strong, while the intermixing of Ti to the Pt phase is 
weaker. Considering the absolute values the predicted 
asymmetry is in excellent agreement with that of given 
by the experiment. On the other hand the experimen- 
tal depth resolution is much larger for Pt/Ti than that 
predicted by TRIDYN. 

It is well known that this code cannot account for 
the relaxation occurring after the collisional cascade [33] . 
Hence we also carried out MD simulations. 

Further advantage of the MD simulations is the ap- 
plied many body tight-binding potentials [34] which are 
known to be more accurate and reliable than the two- 
body potentials used in the TRYDIN code [26,21,33]. 

MD simulations provide a « 8 ML (~ 2.0 nm and 
cr « 16 ML (~ 4.0 nm thick interface after 200 ion im- 
pacts, respectively. The computer animations of the sim- 
ulations together with the plotted broadening values at 
the interface in Fig 3 also reveal the stronger interdiffu- 
sion in Pt/Ti [35]. 

Moreover, the applied setup of the simulation cell, in 
particular the 2.0 nm film thickness is assumed to be 
appropriate for simulating broadening. Our experience 
shows that the variation of the film thickness does not 
affect the final result significantly, except if ultrathin film 
is used (e.g. if less than ~ 1 nm thick film). At around 
5 or less ML thick film surface roughening could affect 
mixing, and vice versa [25]. Also, we do not carry out 
complete layer-by-layer removal as in the experiment. It 
turned out during the simulations that the ions mix the 
interface the most efficiently when they are initialized 
~ 1 ± 0.3 nm above the interface. This value is naturally 
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in the range of the projected range of the ions. Hence, 
the most of the broadening is coming from this regime of 
ion-interface distance . Initializing ions from the surface 
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FIG. 3. The simulated broadening at the interface in nm 
as a function of the number of ions at 500 eV ion energy. 
The ions are initiated from 1.0 nm above the interface. The 
error bars denote a statistical uncertainty in the measure of 
broadening. 

it takes longer and more ions are needed to obtain the 
same level of damage and broadening at the interface as 
it has been found 

by ion-bombarding from inside the film. 

Indeed, we find that the measured ion-sputtering in- 
duced broadening of a « 2 nm for Ti/Pt is in nice agree- 
ment with the simulated value. Hence we expect that 
the most of the measured a is coming from intermixing 
and interface roughening contributes to a only slightly. 
The nice agreement could be due to that the saturation 
of intermixing (broadening) during ion-sputtering is in- 
sensitive to the rate of mixing in the as received sam- 
ples (<7o ~ 1.5 nm interface width including the interface 
roughening in both samples). This is because a sharp in- 
terface and a weakly mixed one (before bombardment in 
the as-received samples) lead to the same magnitude of 
broadening upon ion-sputtering, because the binary sys- 
tems reach the same steady state of saturation under the 
same conditions (ion energy, impact angle, etc.). This 
is rationalized by our finding that during simulations we 
start from a sharp interface (<r ~ 0) and we get a very 
similar magnitude of a than by AES. 

IV. DISCUSSION 

Asymmetric AES depth profiles (when the broadening 
of A/B interface is different from that of B/A) have al- 
ready been observed [16,37]. The asymmetric behavior 
has been explained by the large relative sputtering yield 
of the elements (preferential sputtering) [16,37,42]. In 



the present experiment the relative sputtering yield of 
Ypt/Yn « 0.7 at 500 eV (at 1500 eV we find Y Pt /Y Ti w 
0.9 and a similar rate of broadening) therefore the mech- 
anism is different. The asymmetry of intermixing is also 
known during thin film growth [18,40,41] as well as dur- 
ing concentration gradient driven intermixing processes 
[38]. Moreover, we found high energy ion-beam mixing 
results which report strongly asymmetric interdiffusion 
e.g. in Au/Cu/Au and in Cu/Au/Cu (with thin tracer 
impurity layers sand witched between matrix layers), al- 
though these results have not been discussed in detail 
[39]. However, no reports has been found which present 
asymmetric IM for ion-sputtering with vanishing differ- 
ence of the sputtering yield of the constituents. This is 
interesting result because a new mechanism should be 
worked out to explain the preferential interdiffusion of 
Pt in Pt/Ti. 

The comparison of the measured a with the simulated 
broadening using the 84 — 16 % rule in both cases can 
only be carried out with great care. In principle, these 
values are not comparable directly. However, we make 
some simple assumptions at this point. We expect that 
interface roughening has smaller contribution to broad- 
ening than intermixing. MD simulations supports this 
assumption since we find that the simulated a is in the 
range of the measured values. Unfortunately we have 
no results for the roughening of the samples after ion- 
bombardment. However, we expect that if the simulated 
&mix is comparable with the measured a, than we can 
expect that the ion-sputtered a ro <C <J mix . 

According to the simulations, the mixing of the Pt in- 
creases with the fluence, at 200 ions irradiation (which 
is the highest we reached because of the limited CPU 
time) a w 40 A was found. This is not a saturation value 
(in contrast to the Ti/Pt case) hence we can expect fur- 
ther increase in a. The simulated values of broadening 
are purely coming from intermixing. We expect that the 
rest of measured a is coming from roughening. Neverthe- 
less, MD simulations reproduce the mixing asymmetry 
and we are able to capture the essential features of the 
phenomenon. 

Finally we should also mention that the reason of 
asymmetry is not fully understood yet. We can rule out 
the asymmetry of mass effect with the following computer 
simulation: if we interchange the atomic masses and leav- 
ing all other parameters are unchanged we also get the 
asymmetry of mixing. Hence not mass-anisotropy is re- 
sponsible for the asymmetry of intermixing. It could also 
be that some other parameters, such as e.g. atomic size 
difference could explain asymmetric interdiffusion (there 
is a large difference in atomic volumes between Ti and 
Pt, and also the interaction potential of Ti is strongly 
anharmonic which could cause the observed anomaly of 
mixing). E.g., it does matter if we ion-bombard the Pt 
or the Ti film and the ion-induced injection of Pt to Ti 
is much easier than that of the Ti atoms to the Pt phase. 
This could be due to e.g. the atomic size difference. How- 
ever, the verification of this hypothesis goes beyond the 
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scope of the present paper. 



V. CONCLUSIONS 

Performing AES depth profiling on Pt/Ti and Ti/Pt 
bilayers we have found strong asymmetry of intermix- 
ing depending on the succession of the layers. We could 
reproduce the mixing asymmetry by means of TRIM 
and MD simulations. We get a nice agreement for in- 
terface broadening in Ti/Pt with experiment while for 
Pt/Ti the discrepancy is relatively large. We conclude 
from this that interface roughening might has a signifi- 
cant contribution to broadening in Pt/Ti. In Ti/Pt ion- 
sputtering increases broadening only slightly: in the as 
received sample we find <tq w 1.5 nm, and a « 2.0 nm is 
measured after AES depth profiling. Although, atomistic 
simulations reproduce the main features of interdiffusion, 
the mechanism of the asymmetry of intermixing remains, 
however, unexplained. 
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